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Subject:  Flow  Over  an  Axl symmetric  Body  In  a Cylindrical  Tunnel 
References:  See  page  38. 

Abstract:  The  flow  past  an  axlsynnnetrlc  body  in  a cylindrical  tunnel  is  of 

practical  importance  in  the  study  of  wind  or  water  tunnel  blockage 
effects  and  In  the  analysis  of  high  speed  tube  transportation 
vehicles.  The  present  study  employs  the  Green's  Function  method  to 
analyze  the  flow  past  arbitrary  axlsymmetrlc  bodies  in  cylindrical 
tunnels.  Towards  this  end,  the  appropriate  fundamental  solution 
to  the  Laplace's  equation  is  developed.  This  solution  automatically 
satisfies  the  boundary  conditions  at  the  tunnel  wall  and  at  the  up- 
stream and  downstream  infinities.  The  boundary  condition  at  the  body 
is  formulated  into  an  integral  equation  of  the  second  kind.  A 
numerical  method  is  employed  for  the  solution  of  the  Integral  equation. 
It  is  shown  further  that  the  limiting  case  of  the  flow  in  a tunnel  of 
infinite  radius  agrees  very  well  with  the  previously  known  flow  in  an 
unbounded  medium.  The  present  method  has  been  employed  for  the  analysis 
of  various  body  shapes  of  Interest.  The  results  show  excellent  agree- 
ment in  cases  when  exact  analytical  solutions  are  available.  The 
present  method  has  also  been  found  to  be  computationally  efficient. 
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Nomenclature 


E(m) 

Fo(k) 

Fj^(s,s') 

F2(s,s') 

g(P,P') 

H 

n 

In(x) 

Jm 

Jn(x) 

k 

m 

K(in) 

Kn(x) 

t 

m 

N(P,P’) 

0 

P 

P' 

P„(r,r') 

QnCr.r*) 

r 

8 


tunnel  radius 

Complete  Elliptic  Integral  of  the  second  kind 
finite  Hankel  transform  of  order  zero 
axial  velocity  at  s due  to  a unit  source  at  s' 
radial  velocity  at  s due  to  a unit  source  at  s' 
regular  part  of  the  Neumann  function 
the  Integral  (A-30) 

modified  Bessel  Function  of  the  first  kind  and  order  n 

m-th  positive  root  of  Jj^(w)  ■ 0 

Bessel  Function  of  the  first  kind  and  order  n 

m-th  positive  root  of  J ' (ka)  “ 0 

o 

Complete  Elliptic  Integral  of  the  first  kind 

modified  Bessel  Function  of  the  second  kind  and  order  n 

total  arc  length  of  the  body  along  the  meridian 

modulus  of  the  elliptic  Integrals 

singular  part  of  the  Neumann  function 

order  symbol 

field  point  or  observation  point 
source  point 

function  defined  by  (A-17) 
function  defined  by  (A-28) 
radius 

arc  length  parameter 
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V(P) 

W(P) 

X 

Y„(x) 
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6p(P-) 
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v" 
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r(p,p’) 
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Nomenclature  (cont.) 

single  layer  potential  at  point  P 
double  layer  potential  at  point  P 
distance  |z-z'| 

Bessel  Function  of  the  second  kind  and  order  n 
axial  distance 

angle  between  z-axls  and  the  normal  to  the  body  pointing  Into  the 
fluid 

Dirac  delta  function 
operator  defined  by  eqn.  (21) 
axlsymmetrlc  Laplaclan  operator 
velocity  potential 
Neumann  function 
constants  given  by  (A-18) 
source  strength 
sum  of  series  (31)  and  (32) 


Superscript 

' prime  refers  to  the  source  point,  except  If  used  over  a function 

It  refers  to  differentiation  with  respect  to  the  argument 

Subscript 

r differentiation  with  respect  to  r 


z 


differentiation  with  respect  to  z 
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SECTION  I 
INTRODUCTION 


1.1  Aim 

The  flow  past  an  axisymmetrlc  body  in  a cylindrical  tunnel  Is  of  practical 
Importance  In  the  study  of  wind  or  water  tunnel  blockage  effects  and  In  the 
analysis  of  high  speed  tube  transportation  vehicles.  This  Is  a classical 
problem  for  which  exact  solutions  are  available  only  for  simple  body  shapes 
such  as  spheroids  [1,  2].  Some  of  the  current  approximate  methods  for  the 
solution  of  this  problem  are  applicable  only  to  bodies  of  simple  shape  and 
can  not  deal  with  bodies  with  flat  noses  and  bumps  [6-9].  A few  other 
methods  satisfy  the  boundary  condition  at  the  tunnel  wall  only  approximately 
and  are  not  suitable  for  large  blockages  [11-1^] . The  present  effort  was 
therefore  undertaken  In  order  to  develop  an  efficient  analytical  method 
applicable  to  the  flow  past  arbitrary  axisymmetrlc  bodies  In  cylindrical 
tunnels. 

1.2  Earlier  Studies 

The  earliest  analysis  of  the  flow  over  a body  In  a tube  was  performed  by 
Lamb  [1]  for  a Ranklne  ovoid.  Lamb  obtained  cue  potential  and  stream  function 
for  a source  placed  on  the  axis  of  the  tube.  By  combining  a source  and  a 
sink  of  equal  strength  In  the  presence  of  a uniform  stream,  he  obtained  the 
shape  of  the  ovoid.  Lock  [2]  obtained  an  analogous  solution  for  the  case  of  a 
Ranklne  ovoid  In  a open  Jet  tunnel  of  circular  cross  section.  The  solutions 
of  Lamb  and  Lock  contained  certain  Infinite  series,  which  converged  slowly 
near  the  singularity  and  failed  to  converge  at  the  singularity  Itself.  At 
the  request  of  Lock,  Watson  [3]  derived  alternate  forms  of  the  solution 
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valid  near  the  singularity  with  excellent  convergence  properties.  Subsequently, 

I 

Lock  and  Johansen  [4]  obtained  wind  tunnel  Interference  corrections  for  a 
Ranklne  ovoid  using  Watson's  result. 

The  velocity  potentials  and  stream  functions  for  various  singularities 
Inside  a tube  such  as  ring  and  point  sources,  doublets  and  vortices  were 
derived  by  Levine  [5].  However,  Levine's  results  also  contained  Infinite 
series  similar  to  those  of  Lamb  and  Lock  and,  therefore,  suffered  from  the 

same  convergence  problems.  Alternate  formulations  for  these  potentials  ^ 

Involving  Infinite  Integrals  of  the  modified  Bessel  functions  were  employed 

by  many  authors.  The  Integral  form  for  the  potential  due  to  an  axial  source 

was  used  by  Satlja  [6],  who  obtained  first-order  wall  corrections  for  the 

flow  over  a spheroid;  and  Goodman  [7]  who  applied  slender  body  theory  to  analyze 

the  flow  about  a body  of  revolution  traveling  In  a tube.  Landweber  and  j; 

Gopalakrlshnan  [8]  used  axial  doublets  and  solved  a pair  of  Integral  equations 

for  determining  the  flow  past  a body  In  a tube.  In  a subsequent  work, 

Landweber  [9]  observed  that  the  basic  potential  for  a source  as  used  by 
Satlja  and  Goodman  was  In  the  form  of  a divergent  Integral  and  corrected  It. 

He  also  derived  the  potentials  and  stream  functions  for  axial  and  ring  sources,  i 

t * 

doublets  and  vortices  In  terms  of  Integrals  of.  the  modified  Bessel  functions.  i 

He  used  axial  sources  to  express  the  stream  surface  condition  on  the  body  as  ' 

i 

an  Integral  equation  of  the  first  kind  and  examined  three  methods  of  solution.  j 

It  should  be  noted  that  the  above  methods  use  axial  singularities  In  | 

their  analyses.  It  Is  known  that  methods  using  axial  singularities  can  not  j 

deal  with  bodies  of  complex  geometry  such  as  those  with  flat  noses  and  bumps  [10]. 

The  methods  employing  surface  singularities  are  better  suited  for  these  i 
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classes  of  bodies.  An  obvious  surface  singularity  technique  Is  to  approximate 
the  Infinite  tube  by  a finite  tube  of  sufficient  length  and  distribute  surface 
singularities  on  the  body  as  well  as  the  tube.  The  Klein  and  Mathew  [11] 
method  employs  this  technique  using  ring  vortices,  while  the  Hess  and  Smith 
[12]  method  uses  ring  sources.  The  boundary  condition  at  the  wall  Is  satisfied 
approximately  by  requiring  the  normal  velocity  to  be  zero  at  the  control  point 
for  each  surface  element  on  the  wall  [12].  The  normal  velocity  elsewhere  on 
the  wall  need  not  vanish  and  therefore  the  fluid  appears  to  "leak"  through 
the  wall.  The  effect  of  this  leakage  becomes  significant  as  the  blockage 
ratio  (ratio  of  the  maximum  body  radius  to  the  tunnel  radius)  Increases  and 
the  accuracy  of  these  methods  becomes  poor.  One  technique  to  reduce  this 
leakage  Is  to  require  that  the  Integral  of  the  normal  velocity  at  each  element 
(or,  equivalently,  the  mass  flow  across  each  element)  be  zero  [11,  13]. 

In  all  these  methods,  the  approximation  of  the  Infinite  tube  by  a finite 
tube  requires  some  care.  In  most  cases,  the  tube  Is  continued  a few  body 
lengths  fore  and  aft  of  the  body  so  that  the  velocities  at  the  ends  of  the 
tube  become  sufficiently  uniform.  A recent  analysis  by  Varsomov  and  Halmov  [14] 
employs  source  discs  at  the  ends  of  a finite  tube  to  satisfy  the  uniform 
velocity  condition  at  Infinity.  In  any  case,  these  finite-tube  methods  are 
approximate  and  require  extra  computational  effort  as  well. 

A surface  singularity  method  satisfying  the  wall  condition  exactly  Is 
given  by  Mlloh  [15],  who  used  ring  vortices  and  obtained  the  exact  solution 
for  the  case  of  a prolate  spheroid.  Mlloh' s method  employed  an  Integral 
equation  of  the  first  kind.  Mathew  and  Majhl  [16]  used  ring  vortices  and 
obtained  an  Integral  equation  of  the  second  kind  for  the  boundary  condition  on 
the  body.  Their  method  satisfied  the  boundary  condition  on  the  wall  exactly. 
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1.3  Scope  of  the  Present  Study 

In  the  present  study,  we  employ  the  Green's  Function  method  to  analyze 
the  flow  past  arbitrary  axlsymmetrlc  bodies  In  cylindrical  tunnels.  Towards 
this  end,  we  develop  the  appropriate  fundamental  solution  to  the  Laplace's 
equation.  This  solution  automatically  satisfies  the  boundary  conditions  at 
the  tunnel  wall  and  at  the  upstream  and  downstream  Infinities.  The  boundary 
condition  at  the  body  Is  formulated  Into  an  Integral  equation  of  the  second 
kind.  We  present  a numerical  method  for  the  solution  of  the  Integral  equation. 

This  method  employs  parabolic  surface  elements  with  linearly  varying  source 
strengths.  We  further  show  that  the  limiting  case  of  the  flow  In  a tunnel 
of  Infinite  radius  agrees  very  well  with  the  previously  known  flow  In  an 
unbounded  medium.  Finally,  we  employ  the  present  method  for  the  analysis  of 
various  body  shapes  of  Interest.  The  results  show  excellent  agreement  In 
cases  where  exact  analytical  solutions  are  available. 

In  summary,  the  present  method  satisfies  the  boundary  condition  at  the 
wall  exactly.  It  also  employs  a tunnel  of  Infinite  length  where  the  boundary 
conditions  at  the  upstream  and  downstream  Infinities  are  satisfied  exactly. 

Since  It  uses  a surface  singularity  technique.  It  Is  capable  of  dealing  with 

bodies  of  arbitrary  shape,  unlike  the  methods  employing  axial  singularities 

[6,  8,  9].  Also,  th..  present  method  uses  an  Integral  equation  of  the  second 

kind  and  therefore  Is  better  than  the  methods  using  Integral  equations 

of  the  first  kind  In  several  respects  [9,  15].  Some  of  these  considerations 

are  the  computational  stability,  ease  and  accuracy  of  solution  and  the  existence 

of  the  solution  for  arbitrary  normal  velocities  on  the  body.  The  last  configuration 

Is  of  significance  In  the  analysis  of  the  boundary  layer  on  the  body.  In  all 

the  above  aspects  the  present  method  Is  comparable  to  that  of  Mathew  and  Majhl  [16]. 
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However,  the'.r  formulation  led  to  the  numerical  evaluation  of  double  Integrals 
Involving  considerable  computational  effort.  The  present  method  Is  capable 
of  dealing  with  body  shapes  of  practical  Interest  and  has  been  found  to  be 
computationally  efficient. 
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SECTION  II 
PROBLEM  FORMULATION 


In  this  section  we  present  the  mathematical  formulation  for  the  flow  past 
an  arbitrary  axlsymmetrlc  body  In  a tunnel.  We  review  the  rudiments  of 
potential  theory  In  order  to  obtain  the  solution  of  the  axlsymmetrlc  Laplace 
equation  by  means  of  the  fundamental  solution  appropriate  to  the  region 
Interior  to  the  tunnel.  This  fundamental  solution,  also  known  as  the  Neumann 
function,  satisfies  the  governing  differential  equation  and  the  boundary 
conditions  at  the  wall  as  well  as  the  upstream  and  downstream  Infinities.  The 
remaining  boundary  condition  on  the  body  surface  Is  formulated  Into  an  Integral 
equation  of  the  second  kind.  It  Is  observed  that  the  kernal  of  the  Integral 
equation  has  poor  convergence  near  the  singularity  while  its  convergence  Is 
excellent  elsewhere.  Therefore,  special  formulas,  valid  near  the  singularity, 
are  developed  for  the  kemal  of  the  integral  equation.  These  formulas  are 
rapidly  convergent  and  are  well  suited  for  computation. 

2.1  Potential  Theory 

With  the  advent  of  high  speed  digital  computers  classical  potential 
theory  has  found  wide  use  In  the  solution  of  practical  flow  problems  [10].  This 
theory  expresses  the  solution  of  the  Laplace  equation  for  a given  region  In 
terms  of  certain  singularity  distributions  on  the  boundary  of  the  region.  The 
potential  theoretic  solution  has  two  distinct  advantages:  (1)  it  is  applicable 
to  very  general  shapes  of  the  region  and  (11)  the  resultant  computations  need 
be  performed  only  on  the  boundaries  of  the  region.  This  Is  In  sharp  contrast 
with  the  finite-difference  techniques  which,  even  though  applicable  to  regions 
of  arbitrary  shapes,  require  computations  throughout  the  region.  We  shall  pre- 
sently Introduce  the  concept  of  the  fundamental  solution  to  the  Laplace 
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equation  and  show  how  the  fundamental  solution  and  the  appropriate  surface 
source  distribution  can  be  used  to  solve  the  Neumann  problem. 

Let  be  a closed,  axls5rtmnetrlc  surface  placed  In  an  Infinitely  long 

cylindrical  tunnel,  S2  [Figure  1].  Let  be  the  region  bounded  by 
Consider  the  region  D between  the  body,  and  S2,  the  tunnel.  Let  us  seek 
the  solution  of  the  potential  equation 

V^$(P)  = 0 , P e D , (1) 

where  7^  Is  the  axlsymmetrlc  Laplaclan  operator. 

A fundamental  solution  of  the  above  differential  equation  Is  defined  as 
the  potential  at  any  point  P due  to  a unit  source  at  an  arbitrary  point 
P'  4 P,  such  that  the  boundary  conditions  at  the  tunnel  wall  as  well  as  at 
upstream  and  downstream  Infinities  are  satisfied.  Throughout  this  work  we  will 
follow  the  convention  whereby  the  source  point  Is  denoted  by  a primed  entity 
and  the  observation  point  (or  the  field  point)  Is  denoted  by  an  unprlmed 
entity.  The  only  exception  to  this  convention  would  be  when  the  primed  entity 
Is  a function.  In  which  case  It  denotes  the  derivative  of  the  function  with 
respect  to  Its  argument. 

Suppose  a fundamental  solution  F Is  given,  defined  on  the  region 
DU  = M and  satisfying  the  equation 

7^  r (P,  P')  - -6p(P")  ; P,  P'  E M , (2) 

where  6p(P')  Is  the  Dirac  delta  function. 

Cross-multiplying  F and  $ and  applying  Green's  theorem,  viz: 


< 

(u7^v  - v7^u)dT  - 

3v 

du 

, 

D 3D 


[l.e.,  take  u ■ F,  v ■ In  Equation  (3)],  we  obtain 


I 
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r(P. 


'3n, 


- $(Q0 


.,3r(P,Q') 


3n_^ 


dS,, 


3D 


$(P06p(P")dTp. 


♦(P),  if  P e int.  3D 


(4) 


If  P e ext.  3d 


The  assumptions  we  have  made  In  applying  the  Green's  theorem  are  that  3D 
is  closed  and  smooth  in  the  sense  that  3D  has  a continuously  varying  normal. 
Under  these  conditions,  we  can  assign  the  value  of  %$(P)  for  the  above  in- 
tegral for  a point  P lying  on  the  boundary  3D.  Since  for  P e 3D  the  integrand 
in  the  surface  integral  is  singular,  this  value  of  ^$(P)  for  P e 3D  would 
mean  that  the  surface  Integral  has  a principal  value. 

Thus  Equation  (4)  provides  an  integral  representation  for  the  solution  of 
(1).  Alternatively,  we  could  consider  the  function  defined  by 


♦ (P) 


r(P,  Q')  a(Q')  - 


dSQ. 


(5) 


where  a and  y ai'e  arbitrary  functions,  continuous  on  3D.  By  the  same  argument 
as  before,  we  see  that  the  above  function,  satisfies  Equation  (1),  with  the 
conditions. 


♦ (Q')  - U(Q')  , 

- o(Q')  , Q'  e 3D  . 


(6) 


We,  therefore,  define  the  single  and  double  layer  potential  functions,  respectively  as 
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V(P) 


a(Q')  r(P,  Q")dSQ.  , 


3D 


(7.1) 


W(P) 


Q 


3D 


(7.2) 


It  follows  from  Equations  (6)  and  (7)  that  the  single  layer  potential  V(P) 
Is  useful  In  solving  the  Laplace  equation  when  the  normal  derivative  of  the 
potential  Is  prescribed  on  the  boundary  (Neumann  Problem)  while  the  double 
layer  potential  W(P)  can  be  employed  when  the  potential  Is  prescribed  on  the 
boundary  (Dirichlet  Problem).  It  is  easily  seen  from  equation  (7.1),  that 


V^V(P)  = 0,  for  P e D . 


However,  the  behavior  of  the  potential  on  the  boundary  is  not  evident.  We 
shall,  therefore.  Investigate  V(Q)  for  Q e 3D.  Let  Vj^(Q)  denote  the  limiting 
value  of  V taken  at  a sequence  of  points  in  the  region  interior  to  3D  and 

tending  to  Q e D from  within.  Let  V^(o)  be  the  external  limit  defined  similarly. 
Let  us  indent  the  axis3rmmetrlc  surface,  S,  in  the  neighborhood  of  the  singular 
point,  Q,  by  a semi-toroidal  surface  of  cross  sectional  radius  e.  [Figure  2]. 

Now  the  principal  value  of  the  surface  integral  (7.1)  is  the  limit  as  e -►  0 of 
the  Integral  over  the  Indented  surface.  This  Integral  consists  of  a contribu- 
tion due  to  the  integration  of  oF  over  the  semi-toroidal  surface.  If  Q”"  lies 
on  this  torus,  we  have 
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r(Q.  QO 


7-  • An  e , 


where  r is  the  radius  of  the  body  at  Q.  We  then  have  the  contribution  from 
this  part  to  the  integral  (7.1)  as 


er 


The  above  quantity  tends  to  zero  as  e -*•  0.  A similar  argument  applies  when  the 
external  limit  is  taken.  We,  therefore,  have 


V^(Q)  * V^(Q)  = V(Q) 


(8) 


3V  r 

Repeating  the  limiting  process  for  ^ Lby  indenting  the  surface  towards  the 
Interior  of  3d],  we  find  that  the  contribution  from  over  the  semi- 

torous  is  given  by 


lim 

e 


' . 0421  f 

0 4T7^er  J 


dS 


^ • 2-.“«  - -a-fil  . 


e -►  0 4Tr  er 


Similarly,  when  the  indentation  is  towards  the  exterior  of  3D,  the  above  limit 
is  - ^ , we  therefore  obtain 


3V 


3n, 


OJ22.  + 15L 

3n  ’ 


(9.1) 


3V 


2j;2l  + iL 

2 ^ 3n, 


(9.2) 


We  conclude,  therefore,  that  the  potential  due  to  a source  layer  (single  layer) 
la  continuous  across  the  layer  while  the  normal  derivative  experiences  a Jump 
equal  to  the  source  density.  We  also  observe  that  in  the  Neumann  problem,  the 


interior  limit 


3V 


of  equation  (9.1)  is  specified  and  this  leads  to  the 
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integral  equation  central  to  our  analysis. 


2.2  Integral  Equation 

In  this  section  we  solve  the  Laplace  equation  for  the  region 
between  the  body  and  the  cylindrical  tunnel  by  means  of  a surface  source 
distribution.  We  also  examine  the  boundary  conditions  at  the  various  surfaces 
that  form  the  boundary  of  the  region.  We  note  that  the  boundary  3D  consists  of 
the  surfaces  S^,  and  [Figure  l].  The  surfaces  and  are  at  a 
large  distance  fore  and  aft  of  the  body  and  we  will  take  the  limiting  value  as 
this  distance  tends  to  infinity. 

Let  us  consider  the  potential  due  to  a source  layer: 


V(P) 


J a(QOr(P, 

3D 


(7.1) 


Suppose  the  boundary  condition 


h(Q) 


is  specified  on  3D. 


On  we  have 


h(Q)  - h(s) 


(10) 


where  Q e and  s is  the  arc  length  parameter  on  Sj^.  On  S2,  we  have 


h(Q)  - 0 , Q e S- 


Suppose  we  require 


a(Q)  - 0,  (Q.  Q ) - 0 


(11) 


for  Q e S^,  then  the  boundary  condition  on  is  satisfied. 


-16- 


1 March  1979 
JFrpjk 


On  S,  and  S,  , 
J 4 

We,  therefore 


Taking  a(Q)  = 

and  S, . On  S 
4 


Therefore, 

h(Q) 

Similarly  on  S 


or 

h(Q) 


we  want  the  perturbation  velocity  to  become  zero  as  z ±<*>, 
have 


h(Q)  - 0 , Q e S^, 


(12) 


0 for  Q e S^,  S^,  we  obtain  the  following  constraints  on 
we  have 


h(Q) 


o (Q')  ^ (Q,  Q")dSQ. 


a(Q')  {11m  g (0,  Q')}  dS, 


z -*■ 


Q e . 


(13) 


, we  have 


h(Q) 


3V 


9n, 


a(Q")  {Lt.  U (Q,  0')}  dS-. 


2 -►  00 


Q e 


(14) 
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It  follows,  therefore,  from  Equations  (12),  (13)  and  (14) 


O(Q')  {lira  (o;  Q')}  dS_.  ■=  0 . 

z -►  + ^ 


(14.1) 


The  Neumann  problem  is  solved  if  we  obtain  the  source  distribution  a(Q)  on 
the  body  surface.  This  is  obtained  from  the  boundary  condition  on  the  body. 

Let  a be  the  angle  between  the  z-axis  and  the  normal  to  the  body  pointing 
into  the  fluid  [Figure  2].  Then 


■5 — = - -r—  cosa  + -5—  sina  , 

dn.  9z  dr  ’ 


h(Q)  - h(s) 


From  equation  (9. 1)  , 


-h(s) 


C^(s)  , A ^ ^ 

~2 ' Tz  ^ 3r 


,)  ■ _ .^--1  + 2Tr  cosa  + ^ sinalr^ds'  . 

4 


The  above  integral  equation  should  be  solved  for  the  source  strength 
distribution  0(3).  Here  s is  the  arc  length  parameter  and  I is  the  total 
arc  length  of  the  body  along  the  meridian.  The  source  strength  o(s) 
must  satisfy  the  supplementary  condition  (14.1).  This  supplementary  con- 
dition is  essential  for  the  existence  of  the  solution  to  the  Neumann  Problem 
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and  also  is  of  physical  significance.  We  shall  return  to  this  point  later 
in  the  analysis. 


2.3  Neumann  Function 

It  remains  now  to  find  the  Neumann  function,  r(P,  P')  of  the  above 
analysis.  Since  V(P)  is  the  potential  due  to  a source  layer  (single  layer) 
and  the  region  is  axisymmetric,  r(P,  P')  has  the  physical  interpretation  of  the 
potential  due  to  a source  ring  at  P' , namely  a ring  centered  at  z'  and  of 
radius  r'  kept  inside  an  Infinitely  long  tunnel  of  radius  a. 

Let  us  enumerate  the  conditions  on  r(P,  P") : 


v^rcp,  ?')  = -6p(p')  . 


3r 


0 for  r - a , 


(2) 


3 z 


z i®® 


Since  r(P,  P')  is  a source  of  unit  strength,  the  mass  flow  consideration  gives 
us 

a 

-|^(P,  P')j2TTrdr  - i 2 >z  ' , 

o 

■ y z < z 

We  express  r(P,  P")  as  the  sum  of  two  parts:  the  singular  part  N(P,  P")  and  the 
regular  part  g(P,  P') , l.e.. 


r(P,  P')  - N(P,  P')  + g (P,  P')  , 

with 

7*N(P,  P')  - -6p(P0  , (16) 
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9N  n * 

“ 0 for  r “ a 


(16.1) 


3N  - . 

yi  • 2Trrdr 


(16.2) 


3n 

= 0 for  z -*■  ±«>  , 


(16.3) 


and 


V^g(P,  PO  = 0 , 


l£ 

3r 


= 0 for  r 


(17) 

(17.1) 


2irrdr 


z > z'  , 


(17.2) 


1&  - 


3z  ‘ 2Tfa^ 


-1 

2 


for  z -►  ±«> 


z <z  ' , 


(17.3) 


Let  us  obtain  the  solution  of  Equation  (16).  Expressing  the  Laplaclan 
and  the  Dirac  delta  operators  In  cylindrical  co-ordinates,  we  have 

N + i N + N , - ^ ^ 6(z-z')  . 

rr  r r zz  2TTr 


(18) 


Let  us  define  the  finite  Hankel  transform  117], 


F^(k)  - 


f(r)  J^(kr)rdr  , 


(19) 


for  the  sequence  of  k,  viz.  kj^,  k2,  ...»  given  by  the  roots  of  the  equation 


J '(ka)  - 0 . 

0 


(19.1) 
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The  Inverse  transform  for  (18)  is  given  by  [17], 


f(r) 


2 Z 


We  also  have  the  result  [17], 


a 

' 

(A^f)  J^(kr)r<ir 

0 


a J^(ka) 


f 

r 


r » a 


-k"F^(k) 


where  A 

o 


is  the  operator  given  by 


Let  us  denote  the  transform  of  N by  N,  i.e. 


(20) 


(21) 


N(2;  r',  z';  k) 


N(r,z  ; r',z  ')  Jp(kr)rdr 


(22) 


Now,  applying  the  finite  Hankel  transform  to  equation  (18) , we  obtain 


aJ^(ka)  (N^) 


r “ a 


-k^  N + N - -6(z  - z') 

zz 


J„(kr)rdr  , 


“ -6(z  - z"*) 


Jo(kr') 

2ir 


Invoking  the  Boundary  Condition 


- 0 for  r - a, 


(16.1) 


we  obtain  the  ordinary  differential  equation 


J„(kr') 

S - k^N  - - 6(z  - z')  -V- 


(23) 
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In  order  to  satisfy  the  condition  at  infinity,  viz.  equation  (16.3),  we 
obtain  the  solution  in  the  form, 

z > z" 

z • < z' 

The  requirements  on  N are  that  N be  continuous  and  the  derivative 

^ have  a discontinuity  ^ at  z = z'.  Thus  A,  - A.  and 

2t7  12 


N = A^e 


-k( z - z') 


k(z  - z') 


dz 


z * z. 


dz 


= -2kA, 


z * z 


2tt 


or. 


We  therefore  have 


N(z;  z',  r';  k)  - 


I 


4n 


Air 


I 

! r 


^-k|z  - z'\ 


(24) 


Let  us  now  Invert  N and  the  aid  of  equation  (20)  to  obtain 


N(r,z;r'’,z') 


‘ * Ij  (kr)J  (krO 
^ I ltJo*(ka) 


(25) 
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We  can  now  easily  verify  that  N satisfies  the  Equation  (.612)  and  (16.3). 
Thus 

■1^  2Trrdr  « ^ ^ *Sgn(  z - z')  J 2(ka) 

d z d K O 

0 

The  integral  in  the  above  equation  vanishes  since 

rJ,(kr) 

J^(kr)rdr  = — 

0 


-aJQ'(ka) 


a 

jQ(kr)rdr  * 

0 


Thus,  the  function  N given  by  equation  (25)  satisfies  all  conditions 
(16)  through  (16.3).  Finally  the  solution  of  equation  (17)  is  given  by 


g(P,  P') 


(26) 


as  can  be  easily  verified. 

We,  therefore,  obtain  the  Neumann  function 


r(r,  z;  r', 


_1 


^ I J^(kr)  JQ(krO' 
kJ^2 (ka) 


.(27) 


The  source  distribution  0(Q)  can  now  be  obtained  from  the  integral 
Equation  (15)  and  the  Neumann  function  (27).  We  then  get 


h(8)  - 


- j 0(8)  + 


0(8) 


Fj^(8,  s')  cosa  + F2(8,  s')  slna 


r'ds'  ,(28) 


where  we  have  defined 
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Fj^(s,  s^) 


2Trf  (z.  r; 


r') 


F2(s,  s') 


3r  , > 

2tt  (z,  r;  z , r ) 


It  follows  from  Equation  (27)  that 


F^(s,  s 


')  = - -y  Sgn(2  - z') 


1 + 


J^(kr)  J^Ckr') 


J^^Cka) 


,(29.1) 


F2(s,  s') 


^ j,  e Jj^(kr)  J^(kr') 

a^  k 


(29.2) 


It  should  be  noted  that  the  source  distribution  given  by  equations  (28)  and 
(29)  should  satisfy  the  condition  (14.1).  This  is  an  essential  condition  for 
the  existence  of  the  solution  to  the  Neumann  problem.  It  is  known  from  potential 
theory  that  the  solution  to  the  Neumann  problem  on  a bounded  domain  D exists  only 
when  the  integral,  taken  over  the  bounding  surface,  of  the  values  assigned  to  the 
normal  dei tvative  vanishes  [18].  As  applied  to  fluid  mechanics,  this 
is  the  continuity  condition  (or  conservation  of  mass)  for  the  control 
volume  D.  Rewriting  the  condition  (14.1)  with  the  use  of  equation  (27)  we  get 

I 

o(s')r'ds'  “ 0 , (30) 

0 


where  s'  is  the  arc  length  parameter  on  the  body  and  IL  is  the  total  arc 
length  of  the  body  along  the  meridian.  The  condition  (30)  states  that  the 
areal  sum  of  sources  over  the  body  should  be  zero.  We  note  that  this  con- 
dition is  automatically  satisfied  by  a closed  body.  We  shall  present  here 
a simple  proof  based  on  physical  reasoning:  it  follows  from  our  analysis 
that  half  the  volume  of  fluid  from  the  source  layer  on  the  body  flows  out 
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of  the  body  while  the  other  half  flows  into  the  body.  Since  there  Is  no  flow 

across  the  body  surface  and  the  body  is  closed^  the  net  flow  into  the  body  must 

be  zero.  Otherwise,  the  continuity  condition  for  the  region  inside  the  body 

would  be  violated.  Therefore,  for  a closed  body,  the  areal  sum  of  sources  is 
zero. 

2.4  Summation  of  Series 

The  infinite  series  in  Equations  (29.1)  and  (29.2)  are  rapidly  convergent 
when  the  axial  distance  from  the  source,  |z  - z'|,  is  large.  However, 
convergence  is  slow  near  the  singularity  and  the  series  in  Equation  (29.1)  does 
not  converge  at  all  for  z « z'.  These  defects  can  be  remedied  by  expressing 
these  series  in  terms  of  elliptic  Integrals  and  convergent  series  involving 
1 2 - z'l  and  elementary  functions  of  r and  r'.  This  part  of  the  Investigation 
presents  quite  a few  features  of  mathematical  Interest  and  Is  an  exercise  In 
Itself.  Therefore,  only  the  results  of  the  investigation  are  presented  here  and 
the  reader  can  consult  Appendix  A for  details. 

Let  us  define  the  following  sums: 


„ 2 J^(kr)  Jo(krO 

“1  “ k J^Mka) 


-kx 

„ . Z ® J.(kr)  J„(krO 

“2  k = 

^ J^Mka) 


(31) 


(32) 


Then  for  small  values  of  x > 


. -1 


+ 


0,  the  following  approximations  are  valid: 

a*xE(m)  . ^2* 

ir(l  - m2)  [x2  + (r  + r')^y/^  a 


+ 


”1 
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^8*  ^6 


2304 


384 


480  5040 


(33) 


'2  “ 27rr[x^  + (r  + r')^]^ 


(r^  - - x^)  E(m) 


x2  + (r  - r0‘ 


+ K(m) 


2a 


384 


^^2 

16  ■ 4 


Q2X"  Qox‘‘'' 

32 


48 


aP  18432 


Q4X^  , Qz"**  Qo"* 


768 


384  1440 


(34) 


The  functions  K(m)  and  E(m)  of  the  above  equations  are  the  complete  elliptic 
integrals  of  the  first  and  second  kind  respectively.  The  modulus,  m,  of  these 
elliptic  integrals  is  given  by 


4rr' 


x2  + (r  + r')2  • 


(35) 


The  and  of  equations  (33)  and  (34)  are  functions  of  r and  r' . These 
functions  and  the  coefficients  are  given  in  Appendix  A. 

The  equations  (33)  and  (34)  contain  an  Important  feature  of  our  solution. 
They  enable  us  to  obtain  the  flow  field  in  an  unbounded  medium  (no  tunnel  wall) 
as  the  limit  a of  the  flow  inside  the  tunnel.  Using  equations  (29.1), 
(29.2),  (33)  and  (34)  we  obtain  the  following  limits: 


11m 
a « 


S(m)  (z  - i') 


?^(s,  a")  - - TT(1  - m2)  L(z  - + (r  + r')2JJ/*  ’ (36) 


I 


I 
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lim 

a 00 


F^Cs,  s') 


2Trr[(z  - z')2  + (r  + 


(r^  - r'^  - x^)  E(m) 


(z  - z’')2  + (r  - r')2 


+ K(m) 


The  modulus,  m,  is  now  given  by  Equation  (35)  with  x replaced  by  (z  - z') . 

The  reader  can  verify  that  the  solution  for  an  unbounded  medium  is  indeed  given 
by  the  integral  equation  (28)  along  with  Eqs.  (36)  and  (37),  [19]. 

There  are  various  techniques  for  the  solution  of  the  Fredholm  integral 
equation  of  the  second  kind;  namely,  Eq.  (28),  [18].  One  method  of 
solution  is  to  discretize  the  iategral  equation  using  some  form  of  numerical 

quadrature  scheme,  Simpson's  rule  for  example.  This  technique,  also  known  as 
algebraization,  reduces  Equation  (28)  to  a matrix  equation  in  the  unknowns 
{a(s^),  i = 1,  ...N}  where  N is  the  number  of  intervals  chosen.  The  resulting 
matrix  equation  will  show  strong  diagonal  dominance  which  ensures 
computational  stability.  (This  is  the  same  reason  why  an  integral 
equation  of  the  second  kind  is  preferable  to  that  of  the  first  kind.) 

2.5  Numerical  Example 

The  Integral  equation  (28)  has  been  solved  numerically  by  the  algebraization 
method.  In  this  method,  the  body  surface  is  divided  into  N segments.  The  seg- 
ments are  assumed  to  be  of  parabolic  shape  and  the  source  strength  Is  assumed 
to  vary  linearly  on  each  element.  Simpson's  rule  is  used  for  the  integration 
Indicated  in  equation  (28).  This  procedure  discretizes  Eq.  (28)  into  a matrix 
equation  in  the  unknoims,  {o(s^),  1 ■ 1,  2,  . . . N},  where  a(sj^)  is  the 
source  stiength  at  the  mid-point  of  the  1-th  segment.  Once  the  singularity  dis- 
tribution is  obtained,  the  velocities  and  the  pressure  distribution  on  the 
body  (or  at  any  point  in  the  flow  field)  are  easily  obtained.  The  computational 
method  has  been  prograrmied  in  Fortran  IV  language  for  the  IBM  370/3033 


J. 


(37) 
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processor  at  the  Pennsylvania  State  University. 

The  exact  solution  for  the  Rankine  ovoid  was  employed  to  check  the  cal- 
culations. For  a given  tunnel  radius,  maximum  body  radius  and  the  source 
locations,  the  shape  of  the  ovoid  was  obtained  by  using  Lamb's  results 
[1].  The  poor  convergence  of  the  Infinite  series  given  by  Lamb  was  over- 
come by  using  Watson's  results  [3]  in  the  neighbourhood  of  the  singularities. 
The  resulting  shapes  for  a family  of  ovoids  is  shown  in  Figure  3.  The  present 
surface  singularity  method  was  employed  to  analyze  the  ovoids,  which  were 
divided  into  48  elements.  The  pressure  distribution  obtained  by  the  pre- 
sent method  shows  excellent  agreement  with  the  analytical  solution  [Figure  3] 
in  all  cases.  The  present  method  has  also  been  used  for  various  body  shapes 
of  interest.  These  results  along  with  the  details  of  the  numerical  solution 
will  appear  in  a companion  memorandum. 
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APPENDIX  A 


Sunnnatlon  of  Series  of  Bessel  Functions 


I.  Define 


, -j  • x/a 

n-l  j/aj 


(A-1) 


where  ...  are  the  positive  zeros  of  Jj^(w). 


Consider  the  function. 


[j,<v)Y  (rw/.)  - J (r»/a)Y  (w)]  , 

J-M • 2 ■ »Jo<' 

^ (A-2) 


The  above  function  has  simple  poles  at  J with  residue 

j 

From  equation  (12),  Section  (3.63),  of  Reference  [20], 


Jj^(z)Y2(z)  - J2(z)Y^(z)  - ^ . 


(A-3) 


Sinca  . 0 have  Yj(JJJ^  . (A-4) 

Z m 

Also,  J ’(jj  - 

1 m 2 m 

Substituting  (A-4)  and  (A-5)  Into  (A-3),  It  follows  that  the  terms  of  series 
are  the  residues  at  of  function  (A-2).  By  Cauchy's  theorem,  the  sum  of 

residues  Is 


'1  2Trl 


[jj^(w)YQ(rw/a)  - jQ(rw/a)Y^(w)]  ^ -xw/a 


Jj^(w) 


■J  ’ we 


jQ(r'w/a)dw. 


where  the  contour  of  Integration  consists  of  a semicircle  of  infinite  radius 
lying  on  the  right  half  of  the  w -plane  and  the  Imaginary  axis  indented  to  the 
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right  of  the  origin  by  a small  radius  thereby  excluding  the  simple  pole 
at  the  origin.  It  is  easy  to  verify  that  the  residue  of  the  integrand  at 
the  origin  (*)  is  2.  Therefore,  the  contribution  to  from  the  semi- 

circular Indentation  around  the  origin  is  -1.  Also,  the  contribution  from  the 
Infinite  semicircle  vanishes.  Therefore  the  contribution  from  the 
imaginary  axis  alone  remains  to  be  evaluated. 

Substituting  w = ±lt  for  the  two  halfs  of  the  imaginary  axis,  and 
adding  the  corresponding  integrals  we  obtain 


-1  + - 
ir 


Il(t)Ko(rt/a)  + lQ(rt/a)Kj^(t) 


t sln(f^-)I  (^)  dt  . (A-6) 

a u a 


First  let  us  consider  the  integral 


S ^ I sln(xt/a)t  dt 


(A-7) 


The  above  integral  can  be  expressed  in  terms  of  elliptic  integrals 
by  using  Klrchoff's  formula  [21],  viz. 

ir/Z 


(bt)lQ(ct)  cos  qt  dt 


1 1?^ 


de 


(b+c)^  - 4bc  sin^  6 


(A-8) 


0 '6 
Differentiating  the  above  equation  with  respect  to  the  parameter  q 

(and  interchanging  the  order  of  integration  and  differentation)  we  obtain 

Tr/2 

de 


KQ(bt)lQ(ct)  sin  qt  t dt  - q 


j [q^  + 


(b+c)^  - 4bc  sln^  0]^^^ 


* The  following  results  [ Reference  25]  are  found  useful  in  studying  the  behavior  of 
the  various  Integrands  at  the  origin.  When  V is  fixed,  as  z -M),  we  have 


Jy(*)  - (Hz)  '^r(u+l) 

Yo(z)  ^ I 4n  z 


Y^(z)-  ^ r(v)(iz)’^ 


)ln  z 
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tt/2 


[q^  + (b+c)^] 


3/2 


de 


[l  - m^sin^O] 


3/2  , 


where  m 


4bc 


q^+(b+c)^ 


Taking  b = r/a,  c = r’/a,  q = x/a,  the  integral  becomes, 

2 

ax  E (m) 


'1  r 2 ^ ^ .n2i3/2  2.  ’ 

ttCx  + (r+r')  J (1-m  ) 


with 


m 


4rr' 


x^  + (r+r')^ 


Next  let  us  evaluate  the  Integral 


1 f ^ 

IT  J Ij^(t) 


I^Crt/a)  IpCr't/a)  sin(xt/a)t  dt  . 


The  product  of  the  two  iQ-functions  can  be  replaced  by  the  Sonine- 

Gegenbaur  Integral  [22],  viz. 

ir 

Iq(z)1q(z')  = i J~  Iq(/z^  + z'^  - 2zz'  cosy)  dy 
0 

2 2 2 

Taking  z = r/a,  z'  = r^/a  and  R = r + r’  - 2rr'cosy  and 
Interchanging  the  order  of  integration,  we  get 

,(t) 


J * 

1 I,  f 'iW 

j " J h<'> 

0 “'o 


w 'o<?>  <r>  ‘ ■ 


The  infinite  Integral  with  respect  to  t in  the  above  equation  has 
been  evaluated  by  Watson  [3]  for  values  of  R and  x satisfying  the 
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(A-9) 


(A-10) 


(A-11) 


(A-12) 


(A-13) 


condition 
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< 4a^.  " (A-14) 

Watson's  investigation  is  elaborate  and  consists  of  expanding  the 
product 

IgCRt/a)  sin(xt/a) 

in  ascending  powers  of  t/a  and  integrating  the  resulting  series  term  by 
term.  The  reader  should  consult  Reference  [3]  for  further  details  in  this 
regard.  Using  Watson's  result  and  carrying  out  the  integration  with 
respect  to  Y.  we  obtain 


=1  = -IH- 


a X E(m) 


7r(l-m^)[x^  + (r+rj^)^] 


^2’'  ^2  x^ 

2t  3/2  a ^3  ^ 4 6 ^ 


5 '‘64 
a 


P-x^  4 U„x  P,  P/X^  P-jX^  6 

-1 — + 5 — ) + — ( — § — + -1 * — \ 

24  120  7 '■  2304  384  480  5040  '' 

a 


+ 0 d/a’)  . 

The  coefficients  ^2^  above  equation  are  defined  as 

00 

(t) 


(A-15) 


, .1  r .2„  v£ 

^2n 


U)  ' 


(A-16) 


and  are  evaluated  in  Reference  [3].  The  P^j's  are  functions  of  the  radii 
r and  r^  and  are  defined  as 


TT 


(r^  + r'^  - 2rr'  cosy)”'^^  dy.  (A-17) 


The  functions  P will  be  evaluated  in  section  III. 
n 


L 
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The  integrals  equation  (A-16)  which  are  of  practical  importance, 

were  first  studied  by  Watson  in  Reference  [3].  Watson  Expressed 
these  Integrals  in  terms  of  a rapidly  convergent  infinite  series.  Subsequently 
Smythe  [23]  encountered  these  Integrals  in  his  analysis  of  the  flow  around  a 
sphere  in  a circular  tube  and  evaluated  them  by  a numerical  quadrature 
scheme  using  Weddle's  rule.  Smythe 's  values  for  these  Integrals  were  also 
employed  by  Miloh  [15]  in  his  study  of  the  flow  past  a prolate  spheroid  in  a 
tube.  The  following  values  were  computed  by  this  author  for  the  y-  's  by 


Watson's  method. 


^2  ~ 0.7968241731 
y^  = 1.2004703485 


y^  = 7.4582940625 
y-  - 96.2205196159 

O 


(A-18) 


II.  Consider  the  sum  of  the  series. 


The  terms  of  the  above  series  are  the  residues  at  J of  the  function 

m 


(A-19) 


[jj^(w)Yj^(rw/a)  - Jj^(rw/a)Yj^(w)]  ^ 


-xw/a 


- w Jj.  (rw/a)  e 


(A-20) 


It  is  observed  that  the  residue  of  the  above  function  at  J is 


-Jl(r1m/a)  Y^(JJ  f > J^Crb^  /a)e 


"l’(Jm> 


Now  “ tt/Tj  V this  last  result  into 


I 
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•-  V 


the  above  expression  proves  the  assertion  (A-20) . Therefore  by  Cauchy's 
residue  theorem 

»i 

- if  Ji(w)Yj^(rw/a)  - Jj^(rw/a)  Yj^(w)  ^ -xw/a 

“2  " ■ m I • 2 ‘ 


— oo£ 


where  the  path  of  integration  is  the  imaginary  axis  since  the  integrand 
is  analytic  (well  behaved)  at  the  origin.  Taking  w = ±it  on  the  imaginary  axis 
and  combining  the  integrals  on  the  two  parts,  we  obtain 


1 f [li(t)Kj^(rt/a)  - I^(rt/a)Kj^(t)]  IpCr't/a)  cos(^£) . t dt. 

" J ^ ^-22) 

0 


Consider  the  integral  given  by 


‘^i<?>  ^ 


(A-23) 


The  integral  can  be  evaluated  by  differentiating  Kirchoff's 

formula,  namely  equation  (A-8)  with  respect  to  b.  Thus, 

^/2 


I 


K2(bt)lQ(ct)  cos  qt  t dt 


(b+c) 

[q^  + 


-2c  sin  0 


(b+c)^-4bc  sin^9] 


3/2 


[q^  + (b+c)2]3/2 


[Till 
(b+c)  f 

Jo 


dO 


(l-m^8ln^0)3/2 


■n/2 

“/ 


sin  6d9 
(l-m^sln^0) 


3/2 


8(bc) 


3/2 


[(b+c)  ^ - K(m)}l 

L d-m'^)  (l-m*^)  J 


] 
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Taking  b - r/a,  c - r'/a  and  q - x/a  we  obtain 


2 ,2  2 

r - r'  - X 


L x^  + (r-r')^ 


^ 2Trr[x^+(r+r')^]>S 

where  m Is  given  by  equation  (A-10) • 
Next  we  evaluate  the  Integral 


E(ni)  + K(m) 


i f , 

» J ij(t)  1 


T /Xt.  ^ 


(A-24) 


(A-25) 


By  the  Sonine-Gegenbaur  relation, 

TT 


£)  » i I 

a ’ TT  0 


(Rt/a)  dy  , 


(A-12) 


where 


R^  “ r^  + r'^  - 2rr'  cos  y . 


Differentiating  eqn.  (A-12)  with  respect  to  r, 

• TT 


Ij^(rt/a)  I^Cr't/a) 


f I^(Rt/a) 
ir  R 

J 


(r-r*  cos  y)  dy 


(A-26) 


and  substituting  In  equation  (A-25),  then  Interchanging  the  order  of  Integration,  we  get 

TT 

c-r'  cos  y) 

dy 


4 


2 r (r-r'  cos  y)  f K (t) 

■ ^2  J R J * ^1 


(Rt/a)  co8(— ) t dt. 
a 


0 0 
We  now  expand  the  product  Ij^(rt/a)  co8(^)  t In  ascending  powers  of 
t and  carry  out  the  two  successive  Integrations.  Combining  this  result 
with  the  results  (A-22)  and  (A-24),  we  obtain 
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2Trr[x^  + (r  + 

M/.  Q. 


(r^  -r'^  - x^)E(in) 

[ o :5 + K(n>)] 


2a 


^ < 


- -?  ( 


16 


2 

X + (r-r  ) 


4 ’ ~ 5 ' 384 

a 


V 

32 


v. 

48 


^2*''  ^0*^  9 

+ 1:^7-  - ) + 0(l/a^)  (A-27) 


7 '•  18432  768  384  1440 


The  coefficients  )J2j^  have  been  defined  previously.  The  functions 


Q are  defined  as 
n 


Q„(r,  r’) 


^ r (r  - r'cosY)  (r^  + r'^  - 2rr ’cosy)"''^  dy  • 

Jo 


The  functions  will  be  evaluated  in  section  III. 


III.  In  this  section  the  function  P^(r,r')  will  be  evaluated.  By 


our  definition  : 


P„(r,r')  . - 


TT 

I 


n/2 


(r^  + r'^  - 2rr'  cos  y)  dy  . (A-17) 


We  can  transform  the  above  Integral  to  elliptic  or  pseudo-elliptic  Integrals. 

It  is  easily  seen  that 

r^  + r'^  - 2rr'  cos  y ■ r^  + r'^  - 2rr'(2  cos^  ^ - 1)  ■ (r  + r')^  (1  - h^  cos^  , 


where 


2 4rr' 


2 * 

(rfr')^ 


Hence  equation  (A-17)  transforms  to 

2 


F„(r,r1  (rtr')"  H„  . 


(A-29) 
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where  is  given  by  the  integral 


H - Is 
n 


71 

/ 


-2  -2  1 


n/2 


(1  - h cos  ) dy 


By  change  of  variable,  y “ Tf  - 26,  we  obtain 

7r/2 


H - 

n 

Jr. 


2 2 

(1  - h^  sin  6)  de 


(A-30) 


The  Integrals  obey  the  following  recurrence  relation  [24]; 


H - — (2  - h^)  H , - — (1  - h^)  H , 
n n n-2  n ' n-4 


(A-31) 


This  recurrence  relation  requires  the  values  of  for  tvjo  successive 

even  or  odd  indices  depending  upon  whether  n Is  even  or  odd.  However, 
our  calculations  employ  only  the  even  values  of  n.  Therefore  the  following 
starting  values  would  serve  our  purpose. 

Hq  - w/2  (A-32) 

H2  - TTd-h^  /2)/2 


It  remains  to  evaluate  the  functions  Q . We  have 

n 


Qn  (r.r 


n 


nil 


2 2 

(r-r*  cos  y)  (r  + r'  - 2rr'  cos  y)  dy 


- il± 


/ 

^ t\ 


2 2 2 

(r  + r'  - 2r*  cos  y/2)  (1  - h cos  y/2)  dy  , 


(r-hr*)"'*'^ 


IT 

[ 


2 2 y „Crt.r*y 

(1-h  cos  ^ ) dy 


I 

^ r\ 


cos^  ^ (1-h^  cos^ 
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Substituting  y “ ^ “ 26  we  obtain,  after  some  algebraic  manipulations, 

,-.n+l 

t(r-r')  + (r+r')  (A-33) 

where  is  given  by  equations  (A-31)  and  (A-32). 

It  should  be  noted  that  the  series  and  E2  arose  in  connection 
with  our  work  on  ring  sources  while  the  series  of  Watson  [3]  pertain  to 
point  sources.  It  is  therefore  possible  to  recover  Watson's  results  when 
we  let  the  radius  of  the  source  ring,  namely  r',  become  zero.  We  then 
obtain 


lim 

r'  0 

and 

11m  E„  = a ^^4 
r'  0 ^ 6x 


(A-34) 


where  and  are  the  series  of  Reference  [3].  The  relations  (A-34) 

are  easy  to  verify  and  therefore  their  proof  is  not  presented  here. 
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Figure  3.  ftahklne  Ovoids  in  Tunnel 
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